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The constant density interior Schwarzschild solution for a static, spherically symmetric 
collapsed star has a divergent pressure when its radius R< |Rs = \GM. We show that this 
divergence is integrable, and induces a non-isotropic transverse stress with a finite redshifted 
surface tension on a spherical surface of radius i?o = 3i?^l — For r < Rq the interior 
Schwarzschild solution exhibits negative pressure. When R — Rg, the surface is localized at 
the Schwarzschild radius itself, Rq = Rg, and the solution has constant negative pressure 
p = —p everywhere in the interior r < Rg, thereby describing a gravitational condensate star, 
a fully collapsed non-singular state already inherent in and predicted by classical General 
Relativity. The redshifted surface tension of the condensate star surface is given by Tg = 
A/c/SttG, where Ak = k_|_ — k_ = 2k_|_ = l/Rg is the difference of equal and opposite 
surface gravities between the exterior and interior Schwarzschild solutions. The First Law, 
dM = dEy + Tg dA is recognized as a purely mechanical classical relation at zero temperature 
and zero entropy, describing the volume energy and surface energy change respectively. Since 
there is no event horizon, the Schwarzschild time t of such a non-singular gravitational 
condensate star is a global time, fully consistent with unitary time evolution in quantum 
theory. The p = —p interior acts as a defocusing lens for light passing through the condensate, 
leading to imaging characteristics distinguishable from a classical black hole. A further 
observational test of gravitational condensate stars with a physical surface vs. black holes 
is the discrete surface modes of oscillation which should be detectable by their gravitational 
wave signatures. 
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I. INTRODUCTION 


The endpoint of complete gravitational collapse of a star whose mass exceeds some critical value 
of order of a few solar masses is widely assumed to be a singular state called a ‘black hole.’ This 
hypothesis rests upon the fact that continued gravitational collapse past nuclear densities cannot be 
halted by any conventional equation of state for high density matter. If a trapped surface from which 
light cannot escape is formed, and the energy density of the matter p plus the sum of the principal 
pressures pi is positive 

3 

/3 + ^Pi>0 (1.1) 

i=l 


then the Einstein equations of classical General Relativity are sufficient to prove that gravitational 
collapse must result in a singularity, and black hole formation is inevitable mm- 


The strong energy condition (1.1) holds for all known forms of matter or radiation, but for one 
significant exception, viz. the cosmological vacuum dark energy believed to be responsible for the 
accelerating expansion of the universe. The equation of state of vacuum dark energy is p = —p < 0 
with constant negative pressure, so that p + 2>p = —2p < 0. Indeed, it is just this inequality that results 
in the defocusing rather than the focusing of geodesics, which is the effective repulsive force presumed 
to be causing the expansion of the universe to accelerate [3]. If such an equation of state, violating 


(1.1) were to be realized within the central regions of a star undergoing gravitational collapse, then 
the same effective repulsive force could prevent the formation of a black hole singularity. 

Had cosmological dark energy and the accelerating expansion of the universe not been discovered 
observationally, other considerations lead naturally to this same p = —p < 0 equation of state. In 
classical General Relativity it is the eq. of state of a (positive) cosmological term. In quantum theory 
vacuum fluctuations lead generically to a non-zero cosmological term, resulting in this same eq. of 
state [1]. The vacuum energy density with p = —p due to gluon condensation appears in the bag model 
for hadrons and is necessary to obtain agreement both with QGD sum rules and with experiment [B]. 


Other violations of (1.1) are predicted by quantum theory dependent upon external conditions such as 
in the Casimir effect, now confirmed by a number of laboratory experiments [5]. It is only the extreme 
feebleness of of the gravitational interaction that prevents the gravitational effects of these expected 
violations of the strong energy condition from being directly observed in laboratory conditions. 

Perhaps most persuasive of all, quantum theory applied to black hole spacetimes, as they are usually 
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described, possessing an event horizon and an interior singularity, leads to a number of severe paradoxes 
associated with the validity of unitary time evolution and the preservation of information apparently 
‘lost’ by quantum matter falling into a black hole [7]. Certainly if macroscopic matter disappears into 
a singularity of space and time all predictive power of present physical theories, whether quantum or 
classical is lost as well. Causality and the classical singularity theorems assure us that this common 
fate of both matter and otherwise successful theories such as the Standard Model is inevitable once the 


event horizon is crossed, provided (1.1) holds everywhere within it. Avoiding this fate and removing 
these difficulties consistently within quantum theory and General Relativity requires the abandonment 


of (1.1), and as a result, some modification of both black hole interiors and/or horizons. 

The singularity theorems illustrate the power of arguments and the broad conclusions that may 
be drawn in Einstein’s theory based on simple general bounds or inequalities obeyed by the matter 


stress-energy, such as (1.1), otherwise independently of its detailed composition or eq. of state. When 
spherical symmetry is assumed, it is possible to prove an even more remarkable set of theorems 
about the existence or non-existence of stable fully collapsed stars, relying on apparently even weaker 
conditions on the matter stress-energy. Assuming only isotropic pressures {pi = p for alH = 1, 2, 3), 
a positive energy density that is a non-increasing function of radius, < 0, and matching the 
interior of a star to the vacuum Schwarzschild solution exterior to it, it was shown by Buchdahl [8] 
that the radius of the star must be greater than a certain finite value, R > |Rs, with Rg = 2GM 
the Schwarzschild radius, in order for the interior solution to possess everywhere finite pressure: cf. 
Appendix Subsequently, it was shown that a rigorous lower bound on the central pressure can 
be established under the same assumptions [9]. These theorems do not require the strong energy 


condition (1.1). The existence of these bounds implies that something drastic must happen to the 
interior of a collapsed star even when measurably (i.e. macroscopically) outside of its Schwarzschild 
radius, before a trapped surface forms, at least under the assumptions of spherical symmetry, continuity 
and monotonicity p' <0, entirely within the domain of classical General Relativity. 


Moreover, the Buchdahl and related bounds are established and saturated by comparison to the 
Schwarzschild interior solution with p' = 0, i.e. p = p = constant mm- This constant density interior 
Schwarzschild solution has a divergent central pressure when its radius first reaches the Buchdahl 
bound R = ^Rg- The assumption of strictly constant density, presumed unphysical, together with 
this pressure divergence of the interior Schwarzschild solution has usually been regarded as reasons 
enough to exclude it from further consideration mm- As a result, the behavior of the solution for 
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-R < f-Rs has been little studied or remarked upon in the literature [Hj- However, the existence of the 
Buchdahl and related bounds in which the constant density solution is the limiting case, makes the 
study of the interior Schwarzschild solution relevant to and potentially quite instructive for the general 
case and realistic application to any spherically symmetric fully collapsed self-gravitating mass. 

When one does consider the constant density interior Schwarzschild solution for R < |Rs, some 
of its rather remarkable features quickly become apparent. First, the pressure divergence moves out 
from the origin to a spherical surface of finite radius 

Ro = 3Ryi-^^ >0 for R<^Rs = ^GM (1.2) 


and a new regular solution for 0 < r < i?o opens up behind it, with negative pressure, violating 


(1.1). Second, as the radius of the star approaches the Schwarzschild radius from above R —?■ Rj", this 
region of negative pressure moves outward from the origin and comes to encompass the entire interior, 
since Rq —)• R^ from below. Third, most remarkably of all, in the above limit the entire interior 
becomes one of eonstant negative pressure with precisely the p = —p de Sitter dark energy equation of 
state. Because of the Buchdahl bound and related theorems one can conclude that negative pressure 
and the dark energy ‘quantum vacuum’ equation of state are already inherent in and predicted to 
occur in classical General Relativity quite generally, under sufficiently severe conditions of spherically 
symmetric gravitational compression, and prior to any formation of a trapped surface or event horizon. 

The crucial feature of the pressure singularity at r = of the interior Schwarzschild solution for 
R < |Rs is that the norm of the static Killing vector vanishes at the same point, touching zero in a 


cusp-like behavior at r = Rq, but otherwise remaining positive (c/. Figs.|4]j^ and Figs.[7]|^. In this 
paper we show that the pressure singularity and cusp is in fact integrable through the Komar integral 
formula for the total mass-energy of a stationary configuration |151 fTH] , and results in a distributional 
(5-function in the transverse stress T^g = = p± T\ localized precisely at r = Rq. Thus in 

the limit R —)> R^, Rq —)• Rj the classical interior Schwarzschild solution describes a non-singular 
gravitational condensate star with a physical surface and finite surface tension, proposed in [laiis] 
(see also |19)). in which the thickness of the thin shell quantum surface layer goes to zero. 

The surface tension may be computed in terms of the discontinuity of equal magnitude and opposite 
signed surface gravities k± = at r = Rq, and in fact, this example serves to generalize the 

Israel junction conditions ini ED to a null boundary layer in an unambiguous way. In terms of 


the (redshifted) surface tension and corresponding surface energy to be determined in Sec. IV the 
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differential First Law of Black Hole Mechanics [221 ES] is straightforwardly modified and recognized 
to be a purely mechanical classical relation between volume and surface contributions to the Komar 
energy. There is no temperature or entropy whatsoever associated with the limiting gravitational 
condensate star configuration, as both are identically zero. Moreover, since the Schwarzschild Killing 
time extends throughout the interior and exterior Schwarzschild solution including for Rq = R = Rg, 
and as there is neither a spacetime singularity nor an event horizon, quantum mechanical evolution of 
fields in such a background static spacetime is clearly unitary. 


The paper is organized as follows. In the next section we review the constant density interior 
Schwarzschild solution and show how negative pressures appear when R < ^Rg, summarizing the main 
results of the paper. The succeeding sections starting with Sec. m assemble the necessary formalism 
to analyze the general stationary configuration of matter and gravity through the Komar mass-energy 


integral, energy flux and surface gravity. In Sec. IV we apply the Komar formula to the interior 
Schwarzschild solution and show that its pressure singularity is integrable, and corresponds precisely to 
a positive transverse pressure p±—p>0 and surface tension Tg localized on the discontinuous pressure 
surface. In Sec. |V] we show how the First Law of energy conservation applied to the Schwarzschild 
interior solution may be recognized as a purely mechanical classical relation of volume and areal 
surface energy, the latter determined by the surface tension, at strictly zero temperature and entropy. 


In Sec. two properties of negative pressure gravitational condensate stars with a physical surface 
that may permit them to be distinguished observationally from black holes are discussed, viz. the 
defocusing of null rays passing through the interior leading to different optical imaging properties, and 
the excitation of discrete frequency surface modes detectable by their gravitational wave signatures. 
Sec. Vllj contains our Conclusions and a Discussion of extension of the simple model presented to 
more general situations and its embedding in a more complete theory. 


There are three Appendices. In Appendix the Bnchdahl and related bounds are summarized. 
Appendix contains the mathematical details of how the 5-function distribution in the transverse 
pressure may be obtained by a careful regulation of the pressure singularity by a small parameter 
e in the limit e —)■ O"'', suggesting also how a physical regulator e oc \/h dependent upon quantum 
corrections may enter in a more complete quantum theory. Finally in Appendix the relationship 
between the Israel junction conditions for a spacelike boundary surface as assumed in |18) and the limit 
in which the boundary between the modified de Sitter interior and Schwarzschild exterior becomes 
null is explained. 
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II. INTERIOR SCHWARZSCHILD SOLUTION AND NEGATIVE PRESSURE 


The general static, spherically symmetric line element in Schwarzschild coordinates is 

ds^ = —f{r) dt'^ + (d9^ + sin^ 9 dcj)^) (2.1) 

h[r) ' 

in terms of two metric functions /(r) and h{r). The Schwarzschild time coordinate t is invariantly 
dehned by the existence of a Killing vector such that 


d d 

dt dx>^ 


( 2 . 2 ) 


and the geometry is independent of t. From (2.1) and (2.2) we see that 


- = -gtt = f{r) 


(2.3) 


is a scalar invariant. The radius r is similarly dehned in an invariant geometric manner by the condition 
that A = dvrr^ is the area of the spherical two-surface at hxed r and hxed t. 


The form of the stress-energy tensor of a general static, spherically symmetric distribution of matter 
may be expressed as the diagonal matrix 




^-p 0 0 0 ^ 

0 p 0 0 

0 0 p_L 0 

0 0 0 p_L y 


(2.4) 


in the Schwarzschild coordinates {t,r,9,(p). The three functions p,p and p± are the mass-energy 
density, radial pressure, and tangential pressure respectively. Thus the general, spherically symmetric 
static conhguration of matter and geometry requires hve functions of r in all. 

These hve functions are required to satisfy the Einstein eqs. of classical General Relativity, whose 
independent information is contained in the two components, 


—G^f = \ — h)] = —SttGTK = SirGp (2.5a) 

dr 

G^j. = —K (h — 1) = SttGT’I = SirGp (2.5b) 

rj dr r^ 
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together with the covariant conservation eq. 


^ dr 2f dr r 


( 2 . 6 ) 


which expresses the pressure balance of forces in static equilibrium. It is sometimes convenient to trade 
the two metric functions /(r) and h{r) for the gravitational potential <I>(r) and the Schwarzschild or 
Misner-Sharp mass [H] within a sphere of radius r, m(r), defined by 


/ = , h = l- 


2Gm 


(2.7) 


respectively. From (|2.5a|) —— = so that m(r) is obtained from this by direct integration 

dr 


pr 

m{r) =4:71 r‘^p{r) dr 

Jo 


assuming m(0) = 0. The remaining Einstein eq. (2.5b) and (2.6) thereby become 


, d^ 
h —— = 


Gm 


dr 


dp 

dr 


+ {p + p) 


+ AirGpr 

d^ _ 2(pj_ — p) 
dr r 


( 2 . 8 ) 


(2.9a) 

(2.9b) 


so that <h(r) is the Newtonian gravitational potential in the non-relativistic limit where p,p± <C /? and 
Gm/r"^ <C 1, /i ~ 1 (in units in which c = 1). 


For a star of total mass M and radius R, the metric functions / and h must match the exterior 
Schwarzschild solution in vacuo 


f — u _ 1 2GM _ ^ Rs 

Jextv) — hextv) — 1 — 1 

r r 


R, = 2GM . 


r > R 


( 2 . 10 ) 


where a possible multiplicative constant of integration in /(r) is fixed by the condition that the line 


element (2.1) approach that of flat space Minkowski with the standard interval of time as r —?• oo. In 


addition, for equilibrium the pressure must vanish at r = 72. Thus we have the boundary conditions 


p{R) = 0 , m{R) = M , 


f{R) = HR) = 1 - § 


( 2 . 11 ) 


for the interior solution at the surface of the star. 
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With these conditions the interior solution for r < R is still under determined, and additional 
information about the matter stress tensor must be supplied. Most commonly, perfect isotropic fluid 
behavior is assumed by setting p± = p. Schwarzschild in his second paper [TO] assumed in addition to 
this perfect fluid behavior that the interior has constant density 


P = 


3M 

AttR^ 


( 2 . 12 ) 


and these two additional conditions allow eqs. (2.5)-(2.6) to be solved in closed form. In that case eq. 


(2.8) is integrated immediately to obtain 


dvr 


M 


m{r) = — pr'^ = —^ , h{r) = 1 — , 0 < r < i? 

3 R-^ 


(2.13) 


where we have defined 


= 


SttG _ 2GM Rs 
-^P = 


R^ ■ 


(2.14) 


Then eliminating the gravitational potential function from eqs. (2.9) gives 

dr 


+ Wd (gV + 8.Gpr) = = 0 

dr 2 r 


(2.15) 


where the last equality is valid if p± = p. In view of (2.13), (2.15) may be written in the separable 
form 


dp 

{p + p){8TrGp + H"^) 


r dr 

2(1 - FV2) 


(2.16) 


whose solution is elementary. Integrating from the outer boundary at r = ii where p{R) = 0 to r gives 


p{r) +p _ Vl - H^R^ 
3p{r) + p Vl - 


or solving for the pressure, 


p{r) = p 


Vl - - Vl - 

3Vl - H‘^R^ - Vl - 
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(2.18) 



























and we also have 


p + p = 2p 


Vl - 


3Vl - - Vl - H^r 


(2.19) 


for r < R. Lastly the solution for / or $ obeying the boundary condition (2.11) is easily found to be 


/(r) = i f3\/l - iL 2 i ?2 _ _ ^ 2^2 


> 0 


( 2 . 20 ) 


completing the constant density Schwarzschild interior solution matched to the vacuum Schwarzschild 
exterior solution at r = R. 

Several remarks about this solution bear emphasis. First, its importance is not due to any as¬ 
sumption (at this point at least) that constant density p = p represents a realistic eq. of state for 
high density matter. Rather it represents an extreme situation which can be used as a bound and 
an instructive model for the general spherically symmetric interior solution. Second, since the eqs. 


for p{r) and /(r) are first order, with boundary conditions (2.11), there is no freedom to adjust the 
first derivative f{R)- However this first derivative is also continuous with the exterior Schwarzschild 


solution (2.10), as long as the interior solution (2.18)-(2.20) remains everywhere regular. Third, the 


solution is regular everywhere except for at most one r in the interior where the denominator in (2.18) 


D = 3\/l - H‘^R‘^ - Vl - 


( 2 . 21 ) 


may vanish in the interval r G [0, R]. Fourth, and most importantly, since / = if D = 0 the 

pressure p{r) diverges at the same value as that at which /(r) vanishes. Otherwise /(r) > 0 and the 
interior solution is regular everywhere else with no horizon. 

The solution of D = 0 is given hy r = Rq where 


3 Vl - H^R'^ = Jl-H^R^ 


or Ro = 3RV1-^^ 


( 2 . 22 ) 


which is pure imaginary if R > |Rs. Hence in this case D > 0 and the solution (2.18) remains finite 


everywhere on the real axis [0, R] in the interior of the star. The pressure ( |2.18 ) is everywhere positive 
and monotonically decreasing outward from its maximum at r = 0, and /(r) remains strictly positive 
everywhere in this case. The positive regular functions p{r) and f^{r) (called the redshift factor) are 
plotted in Figs. and for several values of H^R? = R/Rg > | = 1-125. 
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FIG. 1. Pressure (in units of j5) as a function of r (in units of R) of the interior Schwarzschild solution for various 
values of the parameter R/Rs > 9/8 = 1.125. The upper plot shows p{r) for the values R/Rs = 2.50,1.67,1.25 
(brown, orange, red curves) respectively. The lower plot shows p{r) for the values R/Rs = 1.250,1.136,1.126 
(red, green, blue curves) respectively. Note the change of vertical scale in the latter plot (the red curves are 
the same in each) and the rapid increase of the central pressure p(0) as R/Rs approaches the Buchdahl bound 
1.125. 


Now as R 


I Rs from above, the zero of the denominator D (2.22) approaches the real axis at 


Rq = 0 and a divergence of the central pressure p(0) —)• oo appears with /(O) —)• 0 at this same point. 
Hence at the critical value R = ^Rs, the constant density solution (2.18) with p± = p everywhere 
finite strictly ceases to exist. We analyze this divergence of the pressure in the next several sections. 


If nevertheless we consider (2.18)-(2.20) for Rg < R < |i?s the zero of D at Rq given by (2.22) moves 


outward from the origin to finite values of 0 < Rq < R. Then (2.18) shows that a new regular interior 
solution opens up for 0 < r < Rq where D < 0 and p{r) < 0, while /(r) is again positive. 
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FIG. 2. The redshift factor /= as a function of r (in units of R) of the interior Schwarzschild solution for the 
same values of the parameter R/Rg > 9/8 = 1.125 as in Figs. The brown, orange, red, green and blue curves 
are for the values R/Rs = 2.50,1.67,1.25,1.136,1.126 respectively. Note the approach of /s to zero at r = 0 as 
R/Rs approaches the Buchdahl bound 1.125. 


As the star is compressed further and its radius approaches the Schwarzschild radius R —)• Rf from 


outside, (2.22) shows that the radius of the sphere where the pressure diverges and f{Ro) = 0 moves 


from the origin to the outer edge of the star, i.e. Rq —)• Rj, and in that limit the interior solution 
with negative pressure comes to encompass the entire interior region 0 < r < R, excluding only the 
outer boundary at i? = Rg- Finally, and most remarkably of all, since in this limit H^R? = Rg/R —)• 1 


inspection of (2.18) shows that the entire interior solution then has constant negative pressure 


P = -P: 


for 


r < R = Rq = Rg = 2GM 


(2.23) 


with 

/(r) = i (1 - H^r'^) = I h{r) = ^1 - , H = (2.24) 

corresponding to a patch of pure de Sitter space in static coordinates, although one in which gu is | 
its usual value, so that the passage of time in the interior is modified from what would be expected in 
the usual static coordinates of de Sitter space. Typical profiles of the pressure p{r) and /a for values 
of the radius in the range Rg < R < ^Rg and the approach to Rg are shown in Figs. [^and[^ 

The second metric function h{r) for the same values of the radius in the range Rg < R < ^Rg 
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FIG. 3. Pressure (in units of p) as a function of r (in units of R) of the interior Schwarzschild solu¬ 
tion for various values of the parameter R/Rg < 9/8 = 1.125. The upper plot shows p{r) for the val¬ 
ues R/Rs = 1.124,1.087,1.053 (brown, orange, red curves), where the divergence in the pressure occurs at 
Rq/R = 0.106,0.552,0.761 respectively. The lower plot shows p{r) for the values R/Rg = 1.053,1.010,1.001 
(red, green, blue curves), where the divergence in the pressure occurs at Rq/R = 0.761, 0.959, 0.996 respectively. 
For r < Rq the pressure is negative. Note the change of vertical scale in the plots (the red curves are the same 
in each) and the approach of the negative interior pressure p —>■ —p as R approaches the Schwarzschild radius 
Rg from above and Rq approaches Rg from below. 


and the approach to Rg are shown in Fig. When R = Rg = Rq, the exterior region r > Rg is the 
usual vacuum Schwarzschild solution ( 2.10[ ) with an infinitely thin shell discontinuity and jump of the 
pressure at r = Rg, where the horizons / = h = 0 of the interior de Sitter and exterior Schwarzschild 
solutions coincide. Note that both /(r) and h{r) are non-analytic in r at Rq, having a cusp there, and 
being positive on either side oi r = Rq ^ Rg, although r = Rg is a, null surface, there is no interior 
trapped surface, and no causal event horizon. 
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FIG. 4. The redshift factor /= as a function of r (in units of R) of the interior Schwarzschild solution for the 
same values of the parameter R/Rg < 9/8 = 1.125 as in Figs. The brown, orange, red, green and blue curves 
are for the values R/Rg = 1.124,1.087,1.053,1.010,1.001 respectively. Note the approach of the zero of at 
Rq towards R from below as R approaches the Schwarzschild radius Rg from above. 


h 



FIG. 5. The metric function h as a function of r (in units of R) of the interior and exterior Schwarzschild 
solution for the values of the parameter R/Rg = 2.500,1.667,1.250,1.111,1.000 (brown, orange, red, green and 
blue curves) respectively. The minimum of h approaches zero at r = R, as R as R approaches the Schwarzschild 
radius Rg from above. Note the cusp-like behavior in both this and the previous figure. 

Because of the pressure divergence at Rq which first appears at the origin when R = ^Rg the interior 
Schwarzschild negative pressure solution has been little studied or remarked upon in the literature |14j . 
We show in the next two sections by use of the covariant Komar mass-energy integral, containing in 
its integration measure the redshift factor /a which vanishes at the same radius Rq, that the pressure 
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singularity is integrable, and requires relaxation of the p± = p perfect fluid isotropy assumption. 
When i? < the non-isotropic pressure p± ^ p develops a (i-function distribution at r = Rq whose 
coefficient corresponds to a non-zero surface tension, which is therefore inherent in the Schwarzschild 
interior solution. Then the regular negative pressure solution for R < ^Rs and r < iio is recognized 
as a perfectly viable solution, with a surface boundary layer separating the region of negative from 
positive pressures. In the limit R —)• Rf, Rq —)• Rj, it becomes a non-singular classical alternative 
to a black hole for the fully collapsed state. In fact, in this limit the constant density solution found 
by Schwarzschild nearly a century ago becomes essentially the gravitational condensate star solution 
with a surface boundary layer at Rg, proposed in HZHH]. 


III. MASS-ENERGY FLUX AND SURFACE GRAVITY 


Properly interpreting the pressure singularity of the constant density interior Schwarzschild solution 
requires some background and formalism which we provide in these next three sections, starting in 
this section with the Komar mass-energy formula, energy flux and surface gravity. 

In electrodynamics Maxwell’s eqs. = 47r together with Stokes’ theorem allows integration 

of the anti-symmetric field strength tensor on a two-surface S with directed surface element 


[ = dvr [ J'^dS^ = 47r [ J°dU = 47rQ (3.1) 

Js Jv Jv Jv 

thereby expressing the flux integral of the field through the surface in terms of a volume integral for 
the total conserved charge Q. Here dS^ = 5^^dV is a spacelike three-volume element of flat space 
with timelike normal. 

In General Relativity the analogous covariant volume integral over the matter distribution which 
yields the total mass-energy of an isolated stationary system in a general asymptotically flat spacetime 
(not necessarily spherically symmetric) was given first by Tolman [15j and re-derived more geometri¬ 
cally by Komar [16] . The latter derivation relies upon the fact that a stationary spacetime is invariantly 
characterized as one admitting a timelike Killing vector field satisfying the Killing equation 


= 0 


(3.2) 


so that is anti-symmetric in its indices. Thus just as in |3.l|, this anti-symmetric 
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tensor may be integrated over a two-surface and Stokes’ theorem applied to yield 


[ dJ:^= [ R\K^ = AttG [ {2T\-T 5\) d'L^ (3.3) 

J s Jv Jv Jv 

where the commutator of two covariant derivatives in terms of the curvature tensor is used, together 
with ViyK’^ = 0, and the Einstein eqs. = AttG (2T^ — T 5^^ with T = 

It is natural to introduce the time coordinate t with respect to which the spacetime is stationary, 


as in (2.2)-(2.3). If the integration in (3.3) is taken to be a spacelike hypersurface at constant = t, 

(3.4) 


then in curved space 


dV 


where is the vierbein in the orthonormal basis, =diag (—1,1,1,1), and dV is the 

three-volume element in the induced three-metric of the hypersurface. In these natural coordinates 
suited to stationarity, the inverse vierbein = f~ 2 K^ = u^, where is the four-velocity of a 
particle at rest with respect to the time t and / = —K^K^ = —gtt is a spacetime scalar. Thus 

e\ = -r'2K^ = -u^ = fUl (3.5) 

is the directed normal to the hypersurface and = —u^dV = /z 6^^dV. We note that these 
formulae making explicit use of the vierbein frame field and its inverse make sense strictly 
speaking only if is timelike and / > 0, the / —>• 0 limit requiring special care. 


If furthermore the spatial coordinates (x^,x^,x^) of the hypersurface at constant t are defined so 


that the two-surface in (3.3) lies at constant x^, then the directed two-surface element may be written 


dEifiu — e dA 


(3.6) 


with dA the area surface element of the induced two-metric of the surface. The (negative of the) 


integrand of the surface integral in (3.3) is then 


eVe',, = e\) = ^ e^f-'^V^f = E,^ 9^/5 = kj . 


(3.7) 


This quantity is related to the proper four-acceleration of a particle at rest, namely 
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(3.8) 


(jlj^ 1 11 

a}' = — = =-- = — V7 

dr M M j 2f ■’ 


where K'^V^f = dtf = 0 has been used. Multiplying (3.8) by the redshift factor /a converts the 


acceleration with respect proper time r to that with respect to the stationary coordinate time t. 
Taking the projection of lying within the spatial hypersurface in the direction normal to the 


two-surface by contracting with then gives (3.7). Thus k/ is the four-acceleration of the worldline 
of a particle at rest with respect to the time t projected onto the normal to the surface = constant. 


For a three-volume V enclosed by an outer and inner two-surface (9V+ and dV- respectively, (3.3) 
then becomes 

1 


47rG 




(3.9) 


after dividing by dvrG and rearranging. This shows that if the volume V contains no matter, the areal 
surface integral of k/ is independent of the surface dV chosen, and k/ is proportional to the conserved 
mass-energy flux through the surface, analogous to the electric flux normal to the surface of Gauss 


Law (3.1) in electromagnetism. The coefficient l/dvrG has been fixed so that if the surface integral 


over 9V+ is taken outside the matter distribution, it evaluates to the total mass M in the case of an 


asymptotically flat spacetime. Thus in the asymptotically flat case (3.9) becomes 




(3.10) 


expressing the total mass-energy of the system M in terms of a three-volume integral of the matter 
stress-energy, plus a possible surface flux contribution from the inner two-surface. The redshift factor 
/a is then the gravitational redhsift relative to the asymptotically flat region where / = 1. 

If the volume integral over V can be extended to a complete Cauchy surface without an inner 
boundary, such as in the case of a (non-singular) star, then the last surface integral contribution at 
dV- in (3.10) is absent, and ( 3.10| ) gives the total mass of an arbitrary isolated stationary system in 
asymptotically flat spacetime, in terms of a certain volume integral of its stress-energy components. 


The relation (3.10) applied to vacuum solutions of the Einstein eqs. such as the Kerr rotating 


black hole family of solutions, breaks down at the ergosphere boundary where / = 0 and it is no 


longer possible for a particle to remain stationary. In order to extend (3.10) within the ergosphere 


the usual route is to express as a linear combination of the rotational Killing 
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vector = dld(j) and another vector which remains timelike outside the black hole horizon, 

becoming null there. A modified surface gravity can be defined then with respect to this vector, with 
its corresponding redshift factor (—remaining finite and well-defined down to the black hole 
horizon, where it is becomes a constant, kh- The contribution can be evaluated in 

terms of the angular momentum Jh of the black hole. In this way one obtains Smarr’s integral mass 
formula [221E3] 

f^TT 

M = Ah+ 2Q.hJh (3-11) 

for a rotating black hole, with Ah the area of the Kerr black hole horizon and ^h its angular velocity 
of rotation. The differential form of this relation [221123] 


dM = 


Kh 

SttG 


dAn + idn dJn 


(3.12) 


expresses the change of total energy of the system in terms of the change of angular momentum and 
change of surface area of the horizon, and has been called the First Law of Black Hole Mechanics 


The differential form (3.12) suggests that the coefficient 


/ dM \ _ KH 
SttG 


(3.13) 


could perhaps be viewed as the surface tension of the classical black hole horizon j22j . Yet this in¬ 
terpretation is problematic in black hole physics, since if globally extended within its event horizon 
by analytic continuation, it is implicitly assumed that a black hole has no stress-energy whatsoever 
localized on the horizon. Thus it is not clear to what surface energy or surface effect on the horizon 


the ‘surface tension’ (3.13) could possibly be associated. We shall see that the identification of surface 


gravity with surface tension (actually the difference of surface gravities between exterior and interior) 
is possible only when the regular interior solution such as the constant density Schwarzschild inte¬ 
rior solution is known, and when this solution differs fundamentally from that obtained by analytic 
continuation of the exterior vacuum Schwarzschild solution. 


By covariant differentiation of the Killing eq. (3.2), the integrand appearing in (3.10) may be 
written also in the local form 


= 8ttG (rf - 1 5% t'^k^ 


(3.14) 
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prior to any integration. For the general spherically symmetric static form of the metric (2.1) and 


stress-energy tensor (2.4) this eq. takes the form 


A 

dr 


2 jh df 
'V dr 


8ttG r'^ {p + p + 2p±) 


(3.15) 


in Schwarzschild coordinates (2.1). In terms of the surface gravity Hf = k of (3.7), (dropping the 
subscript / henceforth) 


K{r) = 


1 hdl 

2\ f dr 


(3.16) 


so that (3.15) can be written in the form 


^ ^ = 47r (p + p + 2p^) 


(3.17) 


justifying the interpretation of k/G as a mass-energy flux, analogous to the electric flux in Gauss’ 


Law. Indeed since in the exterior vacuum Schwarzschild solution (2.10) 


GM 

Hext[r) — 2 


(3.18) 


the upper limit of the integration of (3.17) outside the matter distribution gives the total mass M. 


IV. SURFACE ENERGY AND SURFACE TENSION 


The local form (3.17) in the spherically symmetric case is the convenient starting point to analyze 


the singular behavior of the pressure (2.18) of the constant density solution at r = Rq. We treat 


the case of general 0 < Ro < R and consider the limit R —)• Rf when Rq —)• Rg at the end. If one 


substitutes the constant density interior Schwarzschild solution assuming p±{r) = p{r) given by (2.18) 
and 

^f = \\D\ (4.1) 


from (2.20) into the right side of (3.17) one finds 
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Att (/9 + 3p) = Airr^ p sgn {D ), r ^ Rq 


(4.2) 


so that the divergence at r = i?o apparently cancels, leaving only a sign function discontinuity 

-1, r <Ro 


sgn {D) = sgn (r - Rq) = 


+ 1, r > Rq 


(4.3) 


at r = Rq. The cancelation of the divergence in (4.2) indicates that the pressure singularity is 


an integrable one with respect to the proper measure in the mass-energy integral (3.10) or (3.17). 
However, since from (2.21) and (4.1), with p oc D~^, and / —)• 0, we are dealing with singular 


distributions rather than smooth functions, and this conclusion is unreliable at the singular point 
itself, potentially missing a local integrable distribution with support only at r = i?o- 


That such a 5-function is indeed present at r = iio is verified by examining the left side of (3.17), 


which is generally valid for all static, spherically symmetric spacetimes. Substituting the interior 


solution (2.13), (2.20) into the quantity to be differentiated on the left side of (3.17) we obtain 


r“ h df 47r _ o 


(4,4) 


The derivative d/dr then produces a 5-function contribution by differentiation of the sign function 
discontinuity, 

A sgn (D) = ^ A sgn (D) = 2 A = 2 5(r - Rq) (4.5) 


(dD\ 

dD 

\ dr / Ro 

dr 


is an even function aX r = Rq. Thus the left side of (3.17) is in toto 


Ro 


l^d^ 

Gdr' 


^r^K(r)^ = dvr sgn(Zl) -|- Rq^ p S{r — Rq) (4.6) 

with a well-defined local 5-function contribution having support only at r = i?o in addition to the 


finite contribution given previously by (4.2). Comparing the general (3.17) with (4.2) and (4.6), the 


5-function distribution must be attributed to the difference 


87r\l ^ {p± -p) = ^pR^^ 5(r - Rq) 


(4.7) 


and hence the breakdown of the isotropic pressure assumption p± = p at the singular radius r = Rq. 
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This interpretation of the 5-function contribution may be confirmed from the pressure balance eq. 


(2.6) expressed in the form 


+ -P)f^ 


(4.8) 


for constant density p = p- Substituting the solution (2.18)-(2.20) in the left side of this relation gives 


dr L‘ 


py/l - sgn {D) = 2pRo Vl - 6{r - R^) = -pR^Jl- ^Rl 5{r - Rq) . (4.9) 

O " 


Multiplying by iirr'^h 2 ^ evaluated at r = i?o then yields (4.7) once again. 


Thus the integrand of the right side of the Komar mass-energy (3.17) may be written in the form 


“V i 


Stt 

p + 3p + 2 {p± - p) = An psgn (D) + — Rq^ p S{r - Rq) 

o 


(4.10) 


in agreement with (4.6). In Appendix [B] we show how the 5-function can be obtained also on the 


right side of (3.15) of (3.17) by a careful regularization of the singularity at r = and taking the 


limit properly of removing the regulator. It is this 5-function contribution, totally integrable within 


the Komar mass formula (3.17), but breaking the assumed isotropic perfect fluid condition p± = p at 


r = Rq, and the fact that the radial pressure is a Principal Part distribution which is also integrable, 
that allows the interior Schwarzschild solution to be interpreted in physical terms, despite failing to 
satisfy the Buchdahl bound tor R < ^Rg. 


From (4.6) or (4.10) we see that the 5-function contribution gives a surface energy contribution 

3 


Eg = ^pR,^ = 2M(^ 


(4.11) 


to the total Komar mass-energy integral (3.10). This is attributable to the discontinuous change of 


sign of the surface gravity (3.7) as Rq is approached from above and below, namely 

GMRq 


4nG 

K± = hm K{r) = ±^— pRo = ± 


R3 


(4.12) 


so that the discontinuity in the surface gravities is 


Ak = k+ — K- = 


2GM Rq RsRq 


R^ 


R^ 


(4.13) 


which leads to a (redshifted) surface tension of the surface at r = Rq of 
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(4.14) 






MRo 


Ak 


2A 8ttRq AttR? SttG 
This is a physicai surface tension associated with a genuine surface energy and positive integrabie 


transverse pressure contribution to the integrai of the Komar mass formuia (3.16). 


The Komar mass formuia may now be consistentiy appiied to the interior and exterior Schwarzschiid 
soiution throughout the domain 0 < r < oo on a compiete Cauchy hypersurface. Since the integrand 


of (3.17) is given by (4.6) or (4.10), the voiume integral of (3.17) excluding r = Rq gives 


r-Rn 


rR 


Ey — — 


I dr + I dr = M — 2M ( ^ 

/o Jr+ V 


(4.15) 


while the J-function at r = i2o gives the surface contribution (4.11) necessary for the total exterior 


Schwarzschiid mass M = Ey + Eg to be obtained. Interestingly, the volume bulk contribution to the 
Komar mass-energy in (4.6) or (4.10) agrees with the Misner-Sharp mass density (2.8) only in 


the outer portion Rq < r < R, where sgn (D) = 1, which is the entire interior region for R > |7?s, 
while for Rg < R < ^Rg and 0 < r < i2o h has the opposite sign. In the latter case the surface energy 


contribution (4.11) is necessary for consistency of the Komar mass formula with the total integrated 


Schwarzschiid mass M from the sum of (4.11) and (4.15). 


V. CONDENSATE STAR LIMIT AND THE FIRST LAW 


In the matching of interior to exterior Schwarzschiid solutions the inner and outer surface gravities 


K± in (4.13) are equal in magnitude, differing only in sign. This is a result of the null surface at r = i2o 


being consistently embedded in a four-geometry, viewed from either side of the surface interface. Indeed 
passing to the limit R —)■ Rf and Rq —)■ Rj , and defining the radial coordinate 




-Rg 


2Rg 


r 


^ Rs 
r 


2 \ 2 


r < Rg 


r > Rg 


(5.1) 


which vanishes at r = Rg, the full interior plus exterior Schwarzschiid solution in the limiting case for 
R = Rq = Rg can be written in the global Rindler-like form 


ds^ = — 


e 


4R? 


dt^ + q\0d^‘^+r^i0dn^ 


(5.2) 
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with 


9(0 = 


r 1 m 


- = f 1 - 

Rl \ mi) ’ 


-Rs < ? < 0 

Q<i<2Rs 


(5.3) 


continuous at r = Rg,^ = 0, and r(^) determined in each region by the inverse of (5.1). Since in this 


form all metric functions are functions of the metric and its first derivative are continuous (C^) 


across the null surface at ^ = 0. In coordinates (5.2) the surface gravity k = Kf defined by (3.7) is 

lei \ 1 


1 d 

K = - 


q \2R, 


2qR, 


sgn (0 


(5.4) 


which because of the continuity of g —)• 1 in (5.3) as ^ 0 is approached from either side, results in 

1 


K± = =h 


2Rg 


(5.5) 


as in (4.12) when R = Rq = Rg. The coordinates (5.2)-(5.3) are admissable in the sense of ref. [25i 


For this matching to the exterior Schwarzschild solution the j factor in the interior de Sit¬ 


ter region (2.20) is essential, which also determines the surface gravities |k+| = |«;_| to being equal 


in magnitude. Thus the limiting case R —>■ Rf of the constant density interior Schwarzschild solu¬ 
tion provides an explicit matching of a (modified) de Sitter interior to the Schwarzschild exterior, 
compatible with general requirements of boundary layers in General Relativity [251 [26] , evading the 
longstanding presumption that such a matching at their mutual Killing horizons is not possible [21] . 
The precise formulation of the matching across a null surface according to the appropriate limit of 
extrinsic curvature tensors is described in more detail in Appendix [C| 


In the limit R = Rq = Rg, (5.2)-(5.5) apply and the discontinuity of the surface gravities and (4.14) 


give the surface tension of the membrane at the null surface r = Rg,^ = 0 


Tg = 


Ak 


K+ — K- 

8ttG 


1 


1 


SirCRg 


(5.6) 


exactly twice the ‘surface tension’ in Smarr’s formula for uncharged, non-rotating black holes (3.11). 


That result takes into account only the outer Schwarzschild geometry while Ak = 2k_|_ is the 


actual surface tension of the transvese pressure term pj_ and surface energy (4.11) associated with 


the pressure jump from interior modified de Sitter to exterior Schwarzschild geometries. The formal 
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FIG. 6. The constant negative pressure p = —p as a function of r/Rg of the limiting form of interior Schwarzschild 
gravitational condensate star solution for R = 1.000001 


Reds hif t 



FIG. 7. The redshift factor as a function of r (in units of R) of the full Schwarzschild gravitational condensate 
star solution for R = Rg. Note the value of the redshift factor /a (0) = 0.5 at the center of the star. 


analogy to a surface tension in Smarr’s black hole formula (3.11) is thus made precise as a mechanical 


property of the physical surface or membrane r = Rg. It is the discontinuity in surface gravities that 


gives a true surface energy contribution (4.11) to the Komar mass-energy from the d-function in p± 


(4.7). Note also that r<j given by (5.6) is small in magnitude for large M. The behavior of the interior 


Schwarzschild solution in the limiting case with the p = —p interior de Sitter solution (modified by 
the I factor in /) as i? —)• R'^ from above and Rq -p- R~ from below is shown in Figs, [ojj^ 

In this (gravitational condensate star or gravastar) limit the volume contribution to the Komar 
mass-energy (|4.15 ) of the p = —p interior is —M and the surface contribution is 2M so that 


E = M = Ey + Eg = -M + 2TgA = -M + 2M . 


(5.7) 
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FIG. 8. The metric function h as a function of r/Rs of the gravitational condensate star solution for R = Rg. 


The surface tension (4.14) is then localized in an infinitely thin shell at the Schwarzschild radius itself, 


and the interior solution is a (modified) static patch of de Sitter space within its horizon, as proposed 
in 


Eq. (5.7) expresses the integral form of the total mass-energy of a gravitational condensate 


star as a volume contribution Ey from the p = —p interior and a surface contribution Eg from the 


thin shell or membrane at r = Rg. It therefore takes the place of the Smarr formula (3.11) for a 
non-rotating black hole. The differential form of (5.7) may also be derived by the methods of [23], by 
which we find 


5M = 5Ey + 2 ts6A+ —^ (5^+ - 5k-) A + (-6M - —^ ) 

8ttG ^ ^ V 47rG J 


(5.8) 


where in the last parentheses the (zero) variation of the Komar energy in the exterior Schwarzschild 
region has been written as a surface integral as in |22j, and (4.14) has been used. Thus the terms 
involving the variations cancel. However since 




(5.9) 


the term with the variation of 5k- adds with 5Ey, and we obtain the differential mass-energy formula 


dM = dEy + Tg dA 


(5.10) 


for the energy conservation (First Law) of a non-rotating gravitational condensate star, in place of 
(3.12) for a black hole. The form of ( |5.10 ) and the previous discussion of the redshifted surface tension 
arising from the transverse pressure p± in the Komar formula now fully justifies the identification of 


the coefficient of 6A in (5.10) as the surface tension of the physical surface located at r = Rg for the 
full (interior -|- exterior) Schwarzshild solution in the limit R ^ Rg. 
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It bears emphasizing that the surface energy contribution Eg in (4.11) and (5.7) does not arise from 
any energy density Tu at the surface, but purely from the transverse anisotropic pressure T ^ = 


p_L which contributes to the Komar energy (3.10). If the Misner-Sharp energy (2.8) is used in (5.7) 
and (5.10) instead, then 6M = p 6Vs = —p dVg with 5Vs = dvr7?^^ SRg, which ascribes all of the energy 
to the volume contribution of constant negative pressure within r < R. Gravitational energy cannot 
be uniquely localized. Only the total energy of the system (matter plus gravitational) is well defined. 
Thus whereas the Misner-Sharp mass m{r) and Komar integrand are quite different locally, their total 
integrals give the same Schwarzschild mass and M = + Eg agree, as seen from the exterior region. 

We close this section with several additional remarks. 

• Since by the Gibbs relation 

p + p = sT + pn (5-11) 

and ^ = 0, no chemical potential corresponding to a conserved quantum number having entered 
our classical considerations, the interior Schwarzschild-de Sitter solution with. p + p = 0 is a zero 
entropy density s = 0 and/or zero temperature macroscopic state, justifying its designation as 
a condensate. 


The differential relation (5.10) expressing the conservation of mass-energy is purely a mechanical 
relation, entirely within the domain of classical General Relativity, rather than a quantum or 
thermodynamic relation. The area A is the geometrical area of the condensate star surface with 
no implication of entropy. There is no entropy at all associated with a macroscopic condensate 


at zero temperature, as (5.11) shows. The Planck length Lpi = ^JhG|(? or Planck mass 
Mpi = sJhcjG or h have not entered our considerations at all up to this point. 

The matching of the metric interior to the exterior solution for R = Rg has the cusp-like 
behavior shown in Figs.[7]j^ which is non-analytic in the original Schwarzschild radial coordinate 
r, invalidating the assumption of complex metric analyticity needed for deriving periodicity in 
imaginary time t. Unlike in the analytically extended vacuum Schwarzschild solution, where 
/(r) becomes negative and the Killing vector becomes spacelike in the interior r < Rg oi the 
black hole, in the negative pressure Schwarzschild-de Sitter interior solution with surface tension 


(4.14) at R = Rg, complex analytic continuation is not possible, and there is no requirement of 
any fixed periodicity in imaginary time. Thus despite its relation to acceleration, the surface 
gravity k carries no implication for temperature or thermal radiation, and its discontinuity is a 
purely mechanical property, namely the surface tension of a physical boundary layer at r = 7?^. 
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Non-analyticity atr = Rs is exactly the property suggested by the analogy of black hole horizons 
to phase boundaries and quantum critical surfaces in condensed matter physics PSIEZ]. Since 
/(r) = 0 corresponds to the ‘freezing’ of local proper time at r = it suggests critical slowing 
down characteristic of a phase transition. The vanishing of the effective speed of light = /(r) 

is analogous to the behavior of the sound speed determined by the low energy excitations at a 
critical surface or phase boundary. This snggests in tnrn that gravitation and spacetime itself 
are ‘emergent’ phenomena of a more fundamental microscopic many-body theory [271128]. 

The positivity of = /(r) on either side of the phase boundary and its interpretation as the 
effective speed of light squared, which must always be non-negative, brings to mind Einstein’s 
original papers on the local Relativity Principle for static gravitational fields, which led him to 
General Relativity from the Minkowski metric = —(? dt^ + -|- dy'^ + dz^, by allowing first 

the time component —gtt = and eventually all other components of the metric to be func¬ 
tions of space (and in general also time) |29|. Thus it could be argued that the non-negativity of 
Cgjj = /(r) in a static geometry, and first order differentiability of the metric in the Rindler-like 


coordinates (5.2), is more faithful to Einstein’s original conception of the Equivalence Princi¬ 


ple, realized by real continuous coordinate transformations, than is complex analytic extension 
around a square root branch point that wonld allow < 0. At the minimnm, the matching 
of the p = —p Schwarzschild interior to exterior provides a consistent logical alternative to 
analytic extension, entirely within the framework of classical General Relativity, provided only 
that surface boundary layers on null boundaries are admitted. 

Finally, since remains timelike for a gravastar, t is a global time and unlike in the analytic 


continuation hypothesis, the spacetime (5.2) is truly static. The t = const, hypersurface is a 


Gauchy surface and is everywhere spacelike. This is exactly the property of a static spacetime 
necessary to apply standard quantum theory, for the quantum vacuum to be defined as the 
lowest energy state of a Hamiltonian bounded from below, and for the Schrodinger equation to 
describe unitary time evolntion, thns avoiding any possibility of an ‘information paradox.’ 

VI. DEFOCUSING OF NULL GEODESICS AND SURFACE OSCILLATIONS 


Since the Schwarzschild time f is a global time, and /(r) does not change sign in the interior 
Schwarzschild solution with R = Rq = Rg, there is no event horizon. The touching of zero of the cnsp 
at r = Rg in Fig. will almost certainly be removed in a more complete theory, such as suggested 
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by the e regulator introduced in Appendix |B| General semi-classical estimates lead one to expect that 
oc h, in which case f{r) would be very small in the vicinity of Rg but nonetheless strictly 


positive everywhere. Even in the limiting case of e —)• 0, where /(r) vanishes at r = Rg, light rays 
with any finite positive radial momentum are still able in principle to pass from the interior outward 
through the Schwarzschild sphere to the exterior region. Since this is very different from the behavior 
of light rays trapped inevitably by a black ‘hole,’ the possibility arises of distinguishing a gravastar 
from a black hole by optical imaging, e.g. by VLBI in the near infrared [301 El]. 

The behavior of light rays in the full p = —p interior plus exterior geometry can be studied by 
means of the geodesic equation 



( 6 . 1 ) 


for zero mass particles = 0. The constants of the motion E and L, energy and angular momentum 
respectively, are defined in terms of the canonical momenta and the affine parameter A along the 
trajectory by 


^ ^ dt 

E^-p, = f(r) - 


(6.2a) 



(6.2b) 


Thus the null geodesic eq. may be written in the form 



(6.3) 


in terms of the effective radial potential 



(6.4) 


which shows that only light rays with vanishing energy E = 0 can hover indehnitely at r = Rg. 


Since angular momentum is conserved, the motion takes place in a plane, which without loss of 
generality can be chosen to be the equatorial plane at 0 = |. Then = L = ^, and by dividing 


(6.3) by for non-radial geodesics, and defining in the usual way the variable 


u = — 
r 


and b = — 
E 


(6.5) 
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the impact parameter, the equation for the null trajectory of a photon may be written as 


du\^ Ih 


- = 


d(j) J 


^ Z u 


AT = 


Rs 


( 6 . 6 ) 


in plane polar coordinates. Substituting the interior solution {2.24) foi R = Rq = Rg we obtain 


du 

d4> 


^Rs 

62 


+ 1 — ^2 , ti > 1 


(6.7) 


du 


whose general solution with — 7 ^ 0 is 

dcp 


4 i ?2 / 4772 

u = \l I + cos{4> - 4>o), r < Rs, I < u < yl + 


in the (modified) de Sitter interior. The upper bound on u corresponds to the radius 


( 6 . 8 ) 


bRs 


^ m.ir). — 


< R.. 


V62 + 4772 

of closest approach to the origin in the interior, which is achieved at (j) = (pQ. At angles 


(6.9) 


4>± = 4>o sin 


the null ray enters and exits the interior region. 


-1 


2R., 


Jb ^ + 4772 


( 6 . 10 ) 


Likewise in the exterior Schwarzschild region one obtains from (6.6) the photon trajectory 

/ du \^ 772 
\d(j)) 52 


= -r^ — +u^ , r > Rs , 0<M<1 


( 6 . 11 ) 


which is solved by the Weierstrass elliptic function V{p)', 53 ) in the form 


u = \+4V{(l);g2,gs) 


( 6 . 12 ) 


in terms of the elliptic invariant parameters 


1 Ri 


92 = 


12 ’ 


93 = 


216 16 62 


(6.13) 
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and we have set a possible second integration constant to zero by choice of the axis at which cj) = 0. 



FIG. 9. The refraction of a null 
ray at the surface of a gravita¬ 
tional condensate star. The an¬ 
gles -dext and -dint are the an¬ 
gles the null ray makes with 
the normal to the surface in 
the Schwarzschild exterior and 
(modified) de Sitter interior re¬ 
spectively. Since according to 
( |6.14D ^ext < the p = -p 
condensate interior behaves as 
a medium with an index of re¬ 
fraction less than unity. 


The important feature of these null geodesics is that they can penetrate and re-emerge from the 

du 


surface at r = By evaluating (6.7) and (6.11) at the boundary u = 1 we see that tani? = 


dcj) 


so 


that they are refracted there according to Snell’s Law 


/ B? 

a/ 1 + ^ 1 + ^ 


(6.14) 


with 'dext and dint the angle the light ray makes with the normal to the surface at r = iig in the 
exterior and interior respectively. This is illustrated in Fig. 

Since dext < dint with equality attained only for radial geodesics with 6 = 0 and dext = dint = 0, 
the condensate interior acts as medium with index of refraction n < 1, or negative lens with respect 
to the vacuum exterior. Hence light rays are defocused by passing through the interior as illustrated 


in Fig. 10, and a gravitational condensate star will have optical imaging characteristics quite distinct 


from a black hole which absorbs all light impinging on its horizon. The detailed imaging expected 
clearly merits a full analysis and modeling in realistic astrophysical environments for comparison to 
observations. 

A second important qualitative difference between gravitational condensate stars and black holes is 
the existence of surface dynamics. In the strictly classical approximation of the interior Schwarzschild 
solution, the surface has zero thickness and additional information about the composition of the surface 
is necessary to determine the normal modes of oscillation. This requires ideally a Lagrangian model 
for the degrees of freedom of the surface, or at the least a phenomenological parameterization of the 
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FIG. 10. The defocusing of null rays passing through the interior of a gravitational condensate star. 


restoring and damping forces acting upon it. A phenomenological treatment of gravitational waves 
due to inspiral was given in ref. [39]. A complete analysis of stability of gravitational condensate 
stars awaits a full dynamical theory. A necessary building block of that more complete theory was 
proposed in |32l |33| , namely the effective action of the quantum conformal anomaly, which possesses 
an additional scalar degree of freedom coupling strongly at the horizon or gravastar surface. 

Even in the absence of a complete theory, it is clear on physical and dimensional grounds that the 
natural frequency of oscillation of the surface modes must be set by the size of the condensate star 
Rs, so that 

ujr. -^ = 101.5 kHz (6.15) 


is the relevant frequency scale. Because the surface is closed, and the surface oscillations will generally 
have a non-zero quadrupole moment, they will generate gravitational waves at characteristic discrete 


quasi-normal mode frequencies of the order of (6.15). This is quite a distinct gravitational wave sig¬ 
nature from the the infalling, coalescence and chirp modes expected for black holes without a physical 
surface [34] , Since the optimal frequency sensitivity of the Advanced LIGO-LSC gravitational wave 
interferometers lie between approximately 50 Hz and 750 Hz, this corresponds to black hole/gravastar 
masses in the range of 10 ^ Mq to a few times 10 ^ Mq 


Further observational tests for distinguishing gravastars with a surface from black holes are possible 
in the case case of rotating gravastars, which have been considered in [36j . and will be taken up in 
more detail future work. 
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VII. CONCLUSIONS AND OUTLOOK 


The constant density interior Schwarzschild solution possesses some remarkable properties that 
seem not to have been appreciated for nearly a century. Its importance lies not in that any ordinary 
matter can be incompressible and exist at constant density under the extreme pressures of gravitational 
collapse, but in the fact that it provides a simple extreme state that saturates the bounds applicable to 
much more general, spherically symmetric configurations; and in the interesting physical insights and 
limiting behaviors it provides. In particular, the Buchdahl bound R > |i?s is based on a comparison 
with the constant density interior Schwarzschild solution. Hence the behavior of the central pressure 
in this solution is expected to be applicable to other equations of state, and to diverge in a similar way 
when a spherically symmetric star with a more realistic eq. of state first contracts to an even larger 
radius outside of its Schwarzschild radius. 

By study of the constant density interior Schwarzschild solution we have shown that negative 
pressure is produced in the center of a spherically symmetric star well before its Schwarzschild radius 
is reached or a trapped surface is formed. Rather than being a reason to reject the solution for R < tRs-, 


we have shown that the pressure divergence is integrable, according to the Komar mass-energy (3.10) 


for a static, spherically symmetric star. A surface energy density localized at the radius (1.2) is 


necessarily produced, with a finite redshifted surface tension given by (5.6) in terms of the discontinuity 
of the surface gravities. The crucial observation is that the static Killing vector becomes null and 
—K^K^ = /(r) vanishes at exactly the same radius as the pressure divergence. Since y/~f multiplies the 


pressure in the Komar energy integral (3.10), the result is that an integrable J-function of transverse 
stress is generated at that radius. With this modification, the isotropy assumption p± = p upon which 
the Buchdahl bound is based is evaded, and the interior Schwarzschild solution becomes again a viable 
model for the non-singular interior of a fully collapsed star, particularly in the limit R ^ Rs- 

As i? —)• Rf from above, the surface discontinuity moves out to the Schwarzschild radius itself, 
i?o j and most remarkably the interior becomes one of uniformly constant negative pressure 

p = —p. Thus the vacuum dark energy eq. of state emerges naturally from classical General Relativity 
under enough spherical compression of matter, and prior to the formation of any trapped surface or 
event horizon. This suggests that although the sequence of constant density conhgurations for R < ^Rg 
may not be physically realistic in detail, the limiting case of R = Rq = Rg may have much broader 
applicability, indicating that a divergence in the central pressure or curvature singularity is avoided by 
a phase transition to a negative pressure condensate and formation of a phase boundary between the 


31 






positive and negative pressure regions. Since in the adiabatic limit, spherical gravitational collapse 
may be conceived as passing through a sequence of slowly decreasing equilibrium states of fixed R, it 
further suggests that the phase transition to p = —p occurs first at the center of the star and moves 
inside out, resulting in an explosive event that would expel prodigious amounts of energy and entropy. 
The final quiescent state of complete gravitational collapse to Rg may be a gravitational condensate 
star with a p = —p modified de Sitter interior, with a hnite surface tension, rather than a black hole. 


Although a fully satisfactory description no doubt requires quantum theory, it is remarkable that 
the possibility, or the prediction, of a phase transition to a negative pressure p = —p equation of state 
counteracting the attractive force of gravity and preventing a singularity in gravitational collapse exists 
already in Einstein’s classical theory, independently of the detailed composition of the matter being 
compressed. This may be less surprising if General Relativity is an effective low energy theory which is 
a limiting case or ‘emergent’ from a more fundamental microscopic description that takes full account 
of the quantum nature of both matter and spacetime |28j . 


From the perspective of the general theory of boundary layers and junction conditions in General 
Relativity, the interior Schwarzschild solution provides an interesting, explicit example showing that 
a second possibility for joining geometries at a null surface, distinct from the analytic continuation 
assumption through a mathematical event horizon usually adopted, is both logically and physically 
possible. Inspection of the line element of the global Schwarzschild solution in Rindler-like coordinates 


(5.2)-(5.3) shows that the spacetime appears to be locally flat (except for the h-fn.) and differentiable 


in the vicinity of the null surface at ^ = 0. The Killing norm 


4R2 


> 0 does not change sign from 


positive to negative as one passes through ^ = 0, and this has the consequence that there is a 5- 
function in transverse stress and curvature localized at ^ = 0. This is very different from the analytic 
continuation hypothesis through complexified {t, r) coordinates, in that Figs. [Tjj^show that the metric 
functions have a non-analytic cusp-like behavior at r = Rg. 


The equality of surface gravity magnitudes |k+| = |k_| (5.5) realized by the explicit interior 
Schwarzschild example shows that the matching of a static de Sitter interior (modified by the \ 
in gtt) to the Schwarzschild vacuum exterior is possible, evading the long presumption that such a 
matching would necessarily produce an unacceptable metric or curvature singularity |21j . Instead 
there is a d-function distributional transverse pressure in an infinitely thin surface layer, leading to a 
finite surface tension and integrable finite surface energy, after accounting for the kinematical redshift 
factor ^y—gtt = \/Y- The non-analytic behavior of the redshift factor, touching zero (classically) at 
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r = Rq, but otherwise everywhere positive, and the sudden change in the vacuum energy from zero 
to p > 0 across the boundary at r = Rg are both strongly suggestive of a quantum phase transition 
and the behavior expected across a quantum critical surface, again pointing towards a more funda¬ 
mental quantum many-body theory of gravitation [28] , and in analogy with similar behaviors in more 
familiar examples in condensed matter systems m The eq. of state of constant p = —p is exactly 
that which should be expected for a quantum macroscopic state with zero entropy and temperature 


from the Gibbs relation (5.11) with zero chemical potential, and hence of a coherent gravitational 
Bose-Einstein condensate of an underlying many-body theory, described here in classical terms. 


In [TH] the factor multiplying the interior de Sitter time in static coordinates was treated as an 
unknown constant C to be determined by the matching through a finite boundary layer. The present 
treatment shows that in order to match properly to the interior de Sitter region, C = | is required in 
the limit of an infinitely thin boundary layer at r = Rg. This results in the interior de Sitter time at 
r = 0 running at half the rate of the asymptotic Schwarzschild time t at infinity. 


The discontinuity of the surface gravities = 2 k+ across the boundary surface also allows 

for a completely mechanical and classical interpretation of the mass-energy changes of a gravastar 


expressed in the First Law (5.10). The surface area A in this relation does not acquire any interpreta¬ 
tion in terms of entropy, and indeed no such interpretation is possible for a condensate interior which 
has zero entropy. Likewise surface gravity acquires no interpretation in terms of temperature, and 
indeed no such interpretation is possible given the non-analytic behavior of the metric functions at the 
cusp. The final macroscopic quantum state attained by a spherical body collapsed to its Schwarzschild 
radius, like that of a neutron star, is one of absolute zero temperature. 


Since the exterior and interior metrics match through (5.2)-(5.3), the Schwarzschild real Killing 
time coordinate t is a global time, and there is no requirement of periodicity in imaginary time which 
would lead to a thermal field theory interpretation, or black hole radiance. The proper interpretation of 


the discontinuity of surface gravities as surface tension in the extended Smarr formula (3.11) and First 


Law (3.12), previously lacking, requires a well-behaved interior solution, and the correct matching 
of the Schwarzschild exterior solution to the (modified) de Sitter interior across a null surface, as 
described in Appendix Since the t = constant hypersurface is a Cauchy surface and everywhere 
spacelike, standard quantum theory can be applied, the quantum vacuum is well defined as the lowest 
energy state of a Hamiltonian bounded from below, and the Schrodinger equation describes unitary 
time evolution, thus avoiding any possibility of an ‘information paradox.’ 
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In this paper we have made no attempt to provide a full theory of the phase transition to the p = —p 
condensate or the surface boundary layer. Instead our aim has been to show how far one can go to 
describe a gravitational condensate, negative pressure interior and surface tension of a fully collapsed 
state completely within Einstein’s classical theory, and therefore consistently with the Equivalence 
Principle, with = /(r) > 0 as Einstein himself first conceived it [29] . 

A more complete quantum treatment will contain as one important element the scalar degree(s) of 
freedom identified in the effective theory of the quantum conformal anomaly |32l [33] , which has been 
shown to be macroscopically relevant, and allow the possibility for the vacuum energy (usually called 
the cosmological ‘constant’) to change [3211371138] . This will almost certainly lead to the infinitesimally 
thin membrane of the classical Schwarzschild solution being replaced by a finite, but still very thin 
surface layer, so that /(r) > 0 strictly. Elementary semi-classical estimates for the physical thickness 
i of that quantum surface layer indicate that I ~ y/LpiRg [IHl |33| . How this physical regulator, 
dependent upon h, might enter is already suggested by the mathematical procedure of regulating the 
(5-fn. distribution in the transverse pressure, discussed in Appendix [B] with e ~ Lpi/Rg. 


Lastly, although difficult to distinguish observationally from black holes, gravitational condensate 
stars do offer some promising possibilities for future astrophysical tests. We have suggested that 
perhaps one of the cleanest tests follows from the defocusing characteristics of the a p = —p interior, 


which in principle can be penetrated by light rays, c/. Fig. If instead of propagating by geometric 
optics through the phase boundary at r = light is scattered in a frequency dependent way, this 
will also produce imaging and lensing characteristics quite different from a black hole. Secondly, the 
existence of a physical surface implies the existence of surface normal modes of oscillation, which 


as a discrete spectrum on the scale of the characteristic frequency (6.15) should be distinguishable 
from the ringdown quasi-normal modes and chirp signals computed assuming there is no surface but 
a mathematical event horizon instead. Clearly the quantitative details of these predictions require a 
fuller theory and more complete treatment, which we defer to a future publication. 


Only when the dynamics of the surface is fully specified can the stability and normal mode spectrum 
of surface excitations be computed reliably. Likewise, only when the interactions of the surface layer 
with ordinary Standard Model matter is fully specified can the question of whether the surface heats up 
or simply absorbs the accreting matter and incorporates it into the condensate interior be addressed. 
The coupling of surface to volume modes and the time scale for the damping of surface oscillations 
are the important quantities needed to obtain quantitative answers to these questions. 
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Appendix A: The Buchdahl and Related Bounds 

A non-trivial bound on the ratio GM/R follows from the classical Einstein eqs. alone, and was 
found by Buchdahl [8], with only the apparently mild assumptions of: 

(i) Isotropic pressure: P± = P 

(ii) Monotonically non-increasing density profile: -j- <0 

dr 

df 

(hi) Continuity of the metric coefficient /(r) and its first derivative — at r = i?. 

Under these assumptions, Buchdahl’s theorem is that miaE] 

r>Irs = Igm (A1) 

o 4 

is the condition for an interior pressure p{r) < oo which is everywhere finite and /(r) >0 non-vanishing 
in the interior r G [0, R\. Conversely, if the three assumptions (i)-(iii) hold, and the star is compressed 
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to a radius R < |-Rs, then the isotropic pressure p must diverge somewhere in the star’s interior. 
In fact, as R ^ divergence appears first at the star’s center, since under the same three 

assumptions (i)-(iii) one can also prove that the central pressure of the star is bounded from below [9] 


p(0) > p 


1 - yi - Rs/R 
3v^l - Rs/R-l 


(A2) 


and the lower bound diverges when R = ^Rg- Moreover the same three assumptions (i)-(iii) are also 
sufficient conditions to prove that the metric coefficient /(r) is bounded from above by [9] 


/(r) < - 3^1 - _ ij 2^2 


1 2 


in terms of the mean density (2.12) with dehned in (2.14). These results show that as R 


(A3) 


iRg 


both the divergence of the central pressure p{0) and the freezing of time /(O) = 0 first occur together 
quite generally at the center of the star, r = 0. 

The inequalities are obtained from the second condition (ii) p' < 0 above, and are saturated by 


the constant density profile p' = 0, p{r) = p, eq. (2.12) for a star of total mass M and radius R. 
Because of these inequalities established for the general spherically symmetric static solution obeying 
conditions (i)-(iii), compared to the constant density Schwarzschild interior solution, the behavior of 
that solution as the Buchdahl limit i? —?■ is approached is of fundamental interest. 

The logical implication of the Buchdahl theorem is that either: 

• The central pressure diverges and a true curvature singularity (black hole) is produced; or 

• One of more of the assumptions (i)-(iii) of the theorem must be relaxed. 


In particular, the way the Buchdahl bound (Al) is evaded by the constant density interior solution 


through the relaxation of condition (i), first at the origin, is instructive and may be expected to be 
quite generic. The study of the behavior of the interior Schwarzschild solution (2.18)-(2.20) for R ^ Rg 


points to the physical resolution, and explicit realization of the second possibility of a non-singular 
interior with an integrable surface energy and surface tension, wherein condition (i) is relaxed. 

Appendix B: Regulating the Singular Distribution at r = Rq 


The naive evaluation of the right side of (3.15) with the singular interior solution (2.18), (2.20) 


gives (4.2), at least for r ^ Rq. The evaluation of the left side (4.6) contains an additional local 6- 
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function contribution. In the Appendix we show that if the singular distributional solution is regulated 
properly, it also gives a (5-function contribution in the limit that the regulator is removed, in agreement 


with (4.6). We treat the case Rq < R first and take the limit R —)■ Rf, Rq —)■ Rg at the end. 


Let us define the dimensionless radial distance x from the singularity of (2.18) by 


r = Ro{l + x). 


(Bl) 


We observe that for Rq < R the singularity in p{r) is a linear divergence in x —)• 0 since the denominator 

(B2) 


(2.21) can be expanded in Taylor series 

D{r) = cix (1 -b C 2 X -I- ...) 


with 


^ dD 
Cl = Rq-j— 
dr 




Ro 




C2 = 


2(1 - H^R,^) 


etc. 


(B3) 


Since the leading term in D ^ is a linear divergence at x = 0, which defines a distribution odd under 
reflection x —)> —x about the singular point, we regularize this divergence by the replacement 


P+P = 


2py/l - H'^R? 2pVl - H‘^R^ 


D 


Cl X 


(l-C2 X-b. . . ^ 


2pVl - H'^R"^ 


X 


Cl 


x^ -b e- 


■^l-C2X-b.. .^ (B4) 


which is the definition of the odd, real Principal Part distribution, V { that is itself integrable. 

\xj 

Likewise for the function <h(r) whose leading singular behavior is even upon reflection x —)■ —x 
about the singular point, we make a similar regulated replacement 


$ = ^ In 
2 


4 


1 


In 


I (x^ -b e^) (l -b 2c2 X -b ...) I = ^ In (x^ -b e^) -b In ^ -b C 2 x -b ... (B5) 

(B6) 


/2 = e'*’ ^ (x^ -b e^) 2 (l + C2 X -b ...) . 


Computing now the derivatives we find 


dp 

dr 


2pVl-H^R^ d 
Rq Cl dx 

2 pVl - H'^R^ 


2^ 2 (1 - C2 a; + • • • ) 

x^ -b ^ 


1 


2x^ 


Rq Cl 


2 C2 X 

X^ -b (x^ -b 6^)2 (x^ -b 6^)2 


-b... 


(B7) 


and 
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yp + p)^ = --(l-C2a^ + ---) ^ ^2 + g2 


dr 


Rq Cl 


^{‘ln(x^ + a+ln(|)+C2X + 


2 pVl - H'^R^ 
Ro Cl 


X 


' / n . o \ o I 


(x^ + (x^ + 6^)2 


(B8) 


in the vicinity of the singular point r = Rq. Summing these contributions and using (2.9b) gives 


2 (px - p) 2pVl - H^R'^ 


Rq Cl (x^ + 6^)2 


(l - C2 X + . . . ) 


(B9) 


so that the leading potentially singular behavior inp_L—p vanishes in the limit the regulator e —?■ 0. This 
explains why it is not seen if the unregulated expressions for p and <h with e = 0 are used uncritically 


as in (4.2). In addition, since the unregulated functions satisfy p± = p, the first subleading term 


proportional to C 2 in (B9) and all subsequent finite terms vanish in the limit e —)• 0 as well. 


Now the combination that appears in the local form of the mass formula (3.15) is 


h 


(fy 

’2(p± -p)' 

[hj 

r 


= AttRq^ p 


[1 - H‘^R^] 5 (1 + xf 

[1 - H^R^ (1 + x)2] 5 (x2 + e2) 


(l + C2 X H-h ) (l - C2 X + . . . ) 


^ ^ 2 - \/l - e2 

= dTri^Q^ p ■ 


y/l - (x2 + e2)i L 


l + 0(x) 


47ri?Q2 p 


(x2 + e2)2 


1 + C>(x) (BIO) 


where (2.13), (2.22), (Bl), (B6), and (B9) have been used. The essential point is that the leading 


term in (BIO) cannot be discarded since although the one parameter sequence of regulated functions 

P2 


1 


= - 


satishes 


2 (x2 + e2)2 


lim (5e(x) =0 for 7 ^ 0 , 

£^•0 


(Bll) 


(B12) 


and thus gives no contribution at any finite x in the limit, it also has the property that 

5^{0) = ^ for any e > 0 


(B13) 


which is just that needed to give rise to a d-function distribution in the limit e —?• 0. Indeed since 

pCSD 

/ dxSe{x) = l Ve (B14) 
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the sequence of as e —)> 0^ defines precisely a (i-distribution 


lim 6e(x) = d(x) 


(B15) 


In view of this result (BIO) becomes 


Svrr^ I'l 


(P± - P) 


YRo^P^(x) 


^R,^p5{r-Ro) 


(B16) 


in the limit e —0, showing how a d-fn. contribution and breaking the isotropic perfect fluid condition 


localized on the singular surface at r = Rq arises on the right side of (3.17), in agreement with (4.6). 


The Komar energy associated with this surface tension is the integral of (B16) or Eg given in 


(4.11). Since pf^ is regular at r = and the integration may be taken over an arbitrarily small 


region surrounding r = Rq, the <5-function contribution and the surface energy in (4.11) is attributable 


entirely to the transverse pressure term p±, and therefore gives rise to a genuine surface tension. 


One subtlety of this derivation is that the Taylor expansion of 74 in (B2) breaks down and ci —)• oo 


at the limiting point Rq = R = Rg, because D defined by (2.21) becomes a simple square root and 


non-analytic in that limit. As a consequence, the metric factor /(r) = vanishes linearly in x at 
r = Rg rather than quadratically as it does for all Rq < Rg. Nevertheless since ci cancels from the 
final combination (BIO), the d-function formula (B16) remains valid in the limit R —)• Rf,Ro —)• R~. 


It is interesting to note that in this limit, even the Principal Part prescription and regulator on the 


radial pressure in (B4) become unnecessary, as p{r) contains no divergence, but only a step function 
discontinuity at r = 7?*, cf. Fig. 


We anticipate also that the formal regulator e introduced in (B4) and (B6) will be replaced by the 


dimensionless radial thickness Ar/Rg of the shell with a small but finite proper width i ~ y/RgAr 

where Ar is expected to be of order Lpi so that e Vh in a more complete quantum treatment. 
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Appendix C: Junction Conditions and the Surface Tension of a Gravitational Condensate Star 

The matching of four-geometries across a general three-dimensional hypersurface interface with 
a J-function distribution of energy density and pressure concentrated on the interfacial boundary 
requires some care, particularly when the boundary becomes lightlike. In this Appendix we show how 
the matching determined by interior Schwarzschild solution for R < |77s, continuity of the metric at 
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r = Ro and the Komar mass-energy may be understood in the context of the general formalism of [40) 
and related to the present authors’ previous work on gravitational condensate stars [TH] • 

In the case when the surface boundary is everywhere timelike, the normal to it n is spacelike and 
may be normalized to unit length 


n • n = n^giMun’^ = +1 


(Cl) 


If the hypersurface of interest is at fixed r in the Schwarzschild coordinates (2.1), this normal is 




(C2) 


The mutually orthogonal basis vectors in the remaining three directions have components 




(C3) 


The extrinsic curvature tensor is then defined by 


Kafe = -rinvR^VuV 


(a) 


rab 


(C4) 


with the Christoffel symbol evaluated in the full four-metric, which is required to be continuous 
(C^) at r = Rq. On the other hand the extrinisic curvature is allowed to be discontinuous at the 
boundary, and its discontinuities 


[K%] = K“+ - K“+ = 47rG 2S% - 6%S' 


(C5) 


determine the surface stress-energy as a 4-function distribution = S\5{r — Rq) concentrated 

on the surface j20j. Indices a,b,c = {t,9,4>} are raised and lowered with respect to the induced 


three-metric on S which, because of (C3) is simply the four-metric of (2.1) restricted to dr = 0. 


Since from (C4) 


^d[_ 

2 / dr 


_ y ■- -j 
^ t— n J- 5 


= K* = ^ 


(C6) 


we have the surface stress-energy components 




47rG 


\/h 


(C7a) 
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S% = s\^a = - 


8ttG 


Vh df 

2 / dr 


+ 


\/h 


(C7b) 


with rj and a the surface energy density and surface tension in the rest frame of the surface at r = Rq- 
From the matching of interior to exterior Schwarzschild solutions or Fig. [^it is clear that there is no 
discontinuity in y^h{r), which would give rise to a surface energy density r/ = —S^f Instead there is a 
rapid dropoff from constant p to zero, which becomes a step function in the limit that the regulator 
of Appendix e —?■ 0 but there is no ^-function in the energy density p, and p = 0. 


From the junction conditions (C6)-(C7) one sees that there is a breakdown of this formalism when 


the surface becomes lightlike due to / —)• 0. Although the discontinuity in ■\Jh{r) may be assumed to 
be small or identically zero so that the surface energy density p vanishes, the surface tension in the rest 
frame of the surface a diverges. This is clearly a simple kinematic effect of the redshift factor /a going 
to zero on a null hypersurface, where an infinite boost would be required and no rest frame exists. In 
the interior Schwarzschild solution this occurs at r = 7?o where / = 0 for any value oi t) < Rq < Rs, 
and in particular persists in the limit Rq —)• Rg, when h{Ro) —)• 0 as well, and the surface coincides 
with the Schwarzschild sphere at r = Rg. 

Because of the apparent difficulty in joining a de Sitter interior to a Schwarzschild exterior at their 
common null boundaries at H = 1/Rg, a maximally stiff surface layer of p = p fluid was assumed to 
be interposed between the de Sitter interior and Schwarzschild exterior solutions in |18j . In that case 
there are two separate timelike boundaries at r = ri and r = r 2 with ri < r 2 both very close to the 


horizon, where the Isreal junction conditions (C6)-(C7) can be applied unmodified. Approached from 
the de Sitter interior side the extrinsic curvature at r = ri in [T^ is 


Vhd 2 2M 1 I 

K/fr, ) =-ln(l — H^r^)\ _ = ——-= 2 =;= —>■ — ■ - - - - ■ — 

* ^ 2 dr ^ y/f - H^rf - H'^rf Rs 


(C8) 


as Hri —)> 1. Since h{ri) = 1 — H^rf is of order e, with Ar = 'r 2 — ri = eRg ~ Lpi this extrinsic 
curvature component and its discontinuity at ri is of order e~ 2 . This accounts for the e~ 2 dependence 
of the surface tensions ai and a 2 at each of the two surface layers in jlHj . However the surface gravity 


K{r) differs from K^*(r) by a factor of f{r) = ^y^h{r) in the Schwarzschild interior solution (2.20), 
so that 


) = t/f{ri ) ^ 




(C9) 
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remains finite, in agreement with (4.12) at Rq = R = Rs- The factor of -\/7(0 is the Tolman redshift 
factor transforming the local extrinsic curvature and surface tension in the surface rest frame to the 
time global Killing time t of the static spacetime. 


Likewise at the outer Schwarzschild boundary we have 


K/(r2+) = —^In (l-^ 
* ^ 2 ^ 2 dr V r 


l^l-Rs/r:^ 


Rs 


(CIO) 


which is again of order e 2 . Multiplying by the exterior Schwarzschild solution redshift factor we 
obtain the surface gravity 

1 


= \/f {ft) = + 


2R, 


= 


(Cll) 


which is again in agreement with (4.12) at Rq = R = Rg. 


Since these surface gravities including the redshift factor tY(r) remain finite in the limit e —>■ 0 as 
the boundary surfaces approach each other at the common de Sitter-Schwarzschild horizon, one can 
take the limits rj” —>■ Rg, dispensing with the interposition of a surface layer of p = p fluid 


entirely and define the surface tension by (4.14), thus obtaining (5.6). In jl8| p oc e 2 —>■ 0 in agreement 


with the results of this work. In this limit therefore one obtains from the interior Schwarzshild solution 
a universal result for the (redshifted) surface tension of a gravitational condensate star, free of any 
model dependent assumptions of the surface layer eq. of state or other properties or parameters. 

This result may be connected with the general formalism of matching across null surface interfaces 
described in m as follows. These authors define a modified or ‘oblique’ extrinsic curvature 


l^ab — 


(C12) 


where N is an oblique or transversal null vector satisfying 


and the limit 


N-N = 0, n-N = —1 


n • n = e —7- 0 . 


(C13) 

(C14) 


is taken. Thus the normal vector n is allowed to approach a null vector as well, in contrast to the 
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uniform spacelike condition (Cl). By a judicious choice of this normalization and the transversal 


vector N and application of the junction conditions to ICab, it is possible to obtain well-defined finite 
results for matching on null surfaces m- Indeed, if one one chooses (for example) 


n'^ = VW K > 


n-n = /~e—)-0 




(C15a) 

(C15b) 


satisfying (C13), one finds 


ICab = -Nf^T^ab = - 


' J Bmfe = / 2 Kafe , 


(C16) 


Note that together with (C6) this implies that any discontinuity of h{r) at the Killing horizon r = Rq, 
where /(r) = 0 would lead to a linear divergence in JC^g = The constant density interior 

Schwarzschild solution avoids this potential problem for any Rq < Rg by the continuity of the metric 
function /i(r) and r] = 0 at r = Rq. 


The modified or oblique extrinsic curvature (C16) has the tt component 


__1 hdl__ 

l^ff - \ I n 1 - ^ 

2^ f dr 


(C17) 


which is discontinuous at r = i? 0 ; leading to the junction condition 


— Tg , A,B — 6,< 


(CIS) 


assuming no discontinuity in ^Jh{r). In the conventions of [iQ] one finds then that 

54 <5(r - Rq) , ^ = rg Hr - Rq) (C19) 

is the distributional stress tensor density on the null surface with only transverse T^g = components. 


The last form of (CI9) is precisely the term in the integrand in the Komar mass-energy (3.15) or (3.17) 


giving a J-function which is integrable with respect to r with the correct measure factor, and (CIS) 


corresponds exactly to the redshifted surface tension of the null surface (4.14) and (5.6), obtained by 


application of the Komar integral formula to the interior Schwarzschild solution. 
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In addition to the redshift factor /a defining the surface gravity and surface tension with respect 


to the Killing time t rather than the proper time of the shell, note that Tg and a of (C7b) have opposite 

signs. Whereas the usual surface tension of a bubble membrane supplies an inward force (corresponding 

to a negative p±) resisting expansion of the enclosed volume with positive pressure, and in the usual 

case the pressure gradient ^ is negative at the boundary, the positive Tg and transverse pressure 

p_L > 0 of the gravitational condensate surface produces an outward force, tending to expand the 

dj) 

surface area, and the pressure gradient is positive at the boundary. Thus we have finally 

dr 


r] = 0 

J.1 

Tg = -f^CT 


(C20a) 

(C20b) 


for the relationship of the present work to ref. [18], with a that obtained by the total discontinuity 
across the full surface layer from Schwarzschild exterior (rj) to modified de Sitter interior (rj~). Using 


(5.4) one may easily check that the same discontinuity (C18) and (5-function distributional transverse 


stress tensor (C19) is obtained in the Rindler-like coordinates (5.2) 


We remark finally that these results, although satisfactory, depend upon the choices (C15) which 


are not unique given the conditions (C13) and (C14). The crucial condition (C15) satisfies is 


n • n = / = -K.Ki^ 


(C21) 


which is the (negative) norm of the timelike Killing vector of stationarity on both sides of the singular 
surface, normalized by referral to the asymptotic time at infinity. It is this normalized Killing field 
that plays a privileged role in the Komar mass-energy and surface gravity 


K = -N^K’' 


-ATMT 




(C22) 


Rescalings n —)• An, N —> A ^N,e —)> A^e, with |A| ^ 1, which would be consistent with ( |C13[ )-( |C14j ) 
in the general formalism of [40] are thereby disallowed by the physical requirement of matching the 
preferred asymptotic stationary Killing time t on both sides of the surface. 


(C13)-(C14 
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